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Abstract 



We consider the AdS/CFT correspondence in the hydrodynamic 
regime up to the second order in a derivative expansion. We demon- 
strate that the fluid conservation equations are equivalent to Ein- 
stein's constraint equations projected on different hyper-surfaces. 
We derive that result for hyper-surfaces of the form r = R (x") up 
to the first order in a derivative expansion of the metric. At the sec- 
ond order expansion, we introduce the notion of different black hole 
horizons, and focus on two particular horizon hyper-surfaces: the 
event horizon and the apparent horizon. We calculate the tempera- 
ture and entropy current for the apparent horizon and show that the 
latter agrees with the area increase theorem for the black hole, and 
differs from the entropy current calculated for the event horizon. 
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Chapter 1 



Introduction 



The AdS/CFT correspondence proposed in [T|, plays a significant 
role in understanding strongly coupled conformal gauge field theo- 
ries. One important application of this duality is the study of the 
effective description of these strongly coupled gauge field theories 
in the long wavelength regini^ This effective description is the 
hydrodynaniic plasma description. The correspondence states that 
a strongly coupled conformal field theory corresponds to a weakly 
coupled string theory, and the latter, in a certain regime, reduces 
to the Einstein classical field equations with a negative cosmological 
constant in an asymptotically Anti de Sitter space-time. This cor- 
respondence is suitable for describing out-of equilibrium fluid flow 
(the hydrodynamic description) which is dual to out-of-equilibrium 
long wavelength dynamics of a black brane. An example for this 
is the fluid-gravity correspondence presented in [2]. The AdS/CFT 
duality has a remarkable application in describing strongly coupled 
QCD type plasma |3]. 

The equations governing relativistic hydrodynamics are conserva- 
tion equations, namely the relativistic navier-stokes equations jUEj. 
Those equations can be found from parts of Einstein's equations 



^We are actually using the small Knudsen number regime, which is explained in the fol- 
lowing chapters. 



in the dual gravitational description [HI [7]. In this thesis we will 
project those equations on different hyper-surfaces and verify that 
every hyper-surface reproduce the same conservation equations. 

From thermodynamics we get the notion of many thermodynamic 
quantities that can be generalized for hydrodynamics, for example: 
temperature, chemical potential and entropy current. The latter will 
be dual to the area of a hyper-surface that admits the area increase 
theorem. In this thesis we will explore the different hyper-surfaces 
and we will focus on two of them, namely: the event and apparent 
horizons. 

We will start by describing the relativistic hydrodynamic fluid 
and the derivative expansion in chapter 2. At the end of the chap- 
ter we will explain the Weyl formalism, which we will work with 
in the rest of the thesis. In chapter 3 we will set the stage for the 
gravity description and in chapter 4 we will provide basic geometric 
quantities. At the end of chapter 4 we will present parts of Ein- 
stein's equation which are dual to the conservation equations of the 
hydrodynamic fluid, and in the following two chapter 5, 6 we will 
explore those equations for different hyper-surfaces. For this we will 
introduce the apparent horizon hyper-surface in chapter 6. In chap- 
ter 7 we will calculate the entropy current and the temperature for 
the apparent horizon, and we will present our conclusion in chapter 
8. 



Chapter 2 



Relativistic Hydrodynamics 



Hydrodynamics j5] is a theoretical model which describes the behav- 
ior of fluids in motion. It treats the fluid as having local domains 
which are in equilibrium, hence one can deflne in those domains 
deflnite thermodynamic quantities such as temperature, pressure, 
entropy, etc. 

One can address a fluid as a relativistic only if the macroscopic 
velocity of the fluid is close to the speed of light. 

2.1 Hydrodynamics equations 

In order to derive the relativistic equations of fluid dynamics we flrst 
need to derive the form of the energy momentum stress tensor. 

If we look on a (i dimensional fluid elemenl|^in its local rest frame, 
we will realize that in that frame Pascal law is valid, which means 
that the pressures on the different surfaces surrounding it are equal. 
Therefore, we get T" = p. On that local frame the component T°° 
of the stress tensor is just the local internal energy density of the 
fluid e. The other components of the stress tensor are zero T°* = 0. 

If we introduce the c?- velocity of the fluid u^, then in the local 
rest frame of the fluid we get u^ = —1 and -u* = . Then the energy 

^Can be understood as a volume within the fluid which contains a large number of molecules 



momentum stress tensor of an ideal fluid will be: 

Tj'l,, = eu^u-'+p{v^'' + u^un (2.1) 

If the fluid equation of state is provided, then one can find the 
connection between the pressure and the energy density. 

From the definition of the d-velocitj in the local rest frame, we 
can calculate its norm (which will be valid in any frame): 

ut'u^ = -1 (2.2) 



We can define the tensor which multiplies the pressure on (2.1) 
as the projector tensor P^'^ , because it projects along the transverse 
direction of the fluid velocity. 

We would like to consider a fluid which undergoes viscosity ef- 
fects, so we will need to introduce a viscous stress tensor t^'^ , which 
will be added to the ideal stress tensor: 

T'' = Klai + r^" (2.3) 

We have some freedom in the determination of the viscous stress 
tensor that we would like to set. First, we will consider the viscous 
stress tensor to be proportional to the derivatives of u^ (x). Second, 
we will need to consider how viscosity will affect the fluid in its local 
rest frame. In order to do so we will require that the energy density 
and the momentum densities in the local rest frame will not change 
due to viscous effects, which means t^^ = 0, r*^* = 0, and because 
the fluid velocity m* = 0, we have in the local rest frame of the fluid 
the following expression: 



r'^^u^ = (2.4) 

This result is correct not only in the rest frame of the fluid but 
also in any Lorentz frame, this result is called the Landau — Lifshitz 



frame. 

We will also specify the equations of motion that can be written 
in a simple form, because they are just the conservation equations 
of the stress tensor: 

d^T^"" = (2.5) 

2.2 Hydrodynamics as an effective Conformal Field The- 
ory 

When working in the regime of small Knudsen number (long wave- 
length) the relativistic field theory has a relativistic hydrodynamic 
description |1]. 

The Knudsen number is just the correlation length of the fluid 
£cor divided by the characteristic length scale L of the variations 
of the macroscopic fields, so we can write the condition of effective 
hydrodynamics description as: 

K^ = ^ < 1 (2.6) 

Under this condition we can expand the stress tensor in a small 
parameter Kn <^ 1 and get: 

oo 

T^"' (x) = Y, Tr {x) , Tr ~ {Kny (2.7) 

1=0 

where Tl^'^ (x) is determined locally by the value of the velocity 
u'^{x), the pressure p(x)p]and their derivatives. For example the 
zeroth order will be the ideal with no derivatives. 

If we consider a relativistic Conformal Field Theory with a finite 
temperature T, we get: 

T; = (2.8) 



^We are assuming that the equation of state has been given. 



Therefore, we can find the equation of state from this condition 
and we get p = ^^ . From dimensional analysis we have p = aT'^ 
and e = (rf — 1) aT"^ where a is a normalization coefficient. Then 



equation (2.1) takes the form: 



Tj'Li = «^' iv'" + du^^n (2.9) 

In order to find the viscous stress tensor t'^", we will use the con- 



formality requirement (2.8), from which we will see that the viscous 
stress tensor has to be traceless, i.e., r'^ = 0. In addition, if we apply 



the equations of motion (2.5) to the zeroth order stress tensor, we 
get a connection between the derivatives of the velocity of the fiuid 
u^ (x) to the derivatives of the fiuid temperature T{x), so we can 
replace, in the next order, the derivatives of T {x) with the deriva- 
tives of M^ (x). We can follow this procedure to all higher order in 
this iterative form, in order to eliminate completely the derivatives 
of T(a;). 

By those two conditions and the symmetry of the stress tensor, 
one can find the stress tensor by taking all possible terms with dif- 
ferent coefficients [8]. Here we will present the complete stress tensor 
to second order in derivatives expansion. 



T^- = aT"" {f]^" + du''u'')-2w'"' + Vru^o) + ^i^i) + ^2^2) + ^3^'"^ 



(2.10) 



where the shear tensor is defined by: 



d u°' 
a,, = P'^^P^d^^up) - -^P,u (2.11) 



d-1 
the vorticity tensor is defined by: 



u,„. = P'^pfjdy^u^^ (2.12) 



'tj.y ^ fi- 



and 



^ ' a — 1 ^ ' a — 1 

Ef2" = 2a>'''cul + 2a''^col Sfg^; = co'^^co"^ - ^^^^^^p^^^. (2.13) 
The brackets represents symmetric or anti-symmetric tensors A(^ai3) 

The coefficients r^, rn, Ai, A2, A3 are called transport coefficients, 
and in the literature [5j 77 is also called the shear viscosity coefficient. 



2.3 Weyl covariant formulation 

In this section, we will introduce the Weyl covariant formulation as 
done in [9J, which is useful for conformal fluid. 

We will consider a c? > 3 conformal fluid in a curved background 
with metric S'^^F] and we will look at some observables of the fluid. 
We will consider a conformal transformation: 

9,u = e'%,, g^'^ = e-'^r'' (2.14) 

where (f) = (f) (x'^). 

As in the previous section we will denote by u'^ the (i-velocity 



of the fluid with the normalization (2.2). From this normalization 



and (2.14) we can flnd the transformation rule of the ^-velocity by: 
g^yU^u'^ = 'g^yU^u^ = —1 and we get: u^ = e~'^u'^. 

However, if we look on the transformation rule of the covariant 
derivative of ■u'^ we get: 



•^In the rest of the thesis we will consider a flat background metric r;,,^, thus for a flat 
background one needs to replace 3^1/ — ^ rj^i/ and V^ — > 9a- 



V^u" = e-^ (V^u" + 6';Wd,<P - 5^^aSV9^0) (2.15) 

which does not transform homogeneously. In order to deal with ho- 
mogeneous transformations, which are allowed in a conformal fluid, 
we define a Weyl covariant derivative Dx in the following manner: if 
a tensorial quantity Q'^^Z transforms Q'^^Z = e~^'^Q'^'.'.'. then the weyl 
covariant derivative will transform:DAQl^.'.'.' = g~^'^DxQ'^Z, where: 



DxQ>:::. = VxQ'iz + coAxQ'^::. 

+ [gxaA^ - 6^,A^ - 6^Ax] Qf:. +■■■ (2.16) 
- [gx^a^ - 51A, - S^] Q^l - . . . 

Where Ax is a one-form that can be determined uniquely by the 
requirements that the Weyl covariant derivative of the rf-velocity 
will be traceless and transverse to the fluid direction, i.e., D^u^ = 
and u^DxU^ = 0. We get: 

A, = u^VxUu - -^fZl""- (2.17) 

Note that the Weyl covariant derivative is metric compatible, i.e: 

Dxg^ly = 0. 

We will present here additional transformations of some observ- 
ables of the fluid: 



= < + < = e-^5^«^ (2.18) 

= - {D^'u" + D^u^") = e-^'^a^"', (2.19) 

- {D^'u" - D^u^") = e-^'^u^'' (2.20) 



w^" 



^A< = (Va - (rf - 1) Aa) a\ (2.21) 



DAC^^^ = (VA-(rf-3)AA)a;^^ (2.22) 



u^Dy^a^, = P^^PfVxa^fs + ^^^^. (2.23) 

We will also define a Weyl covariant Riemann curvature tensor 
for a vector field with a weight w ,i-e-^'' = e~'^'^V'^ by 



[D^, D,] Vx = cuJ^^uVx + n^,^V^ (2.24) 

with 

J-^, = V^A, - V,A^ (2.25) 



T^fiuXa = R/iuXa + J^/iugXa — 4:6['^g^][x6^-^ ( V^A/j + A^A^ 



A^A- 



-QaP 



2 

(2.26) 

After defining the Weyl covariant Riemann curvature tensor, it 



is natural to define the Weyl-covariantized Ricci tensor TZ^i, and 
Weyl-covariantized Ricci scalar TZ by, 



(2.27) 



n = R + 2{d-l) VxA^ -{d-2){d- 1) AxA^ (2.28) 

Note that we will work with a fiat metric on the boundary, hence, 
the Riemann curvature tensor, the Ricci tensor, and the Ricci scalar 
of the boundary are zero (However, their Weyl definitions are not 
zero). 

In order to evaluate fluid dynamics to the second order, we will 
need to take all sorts of terms with two Weyl covariant derivatives. 
This work was done in p!0] and can also be found in [9] which we 
use for our notations. 
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Chapter 3 



Fluid Gravity Correspondence 



The fluid gravity correspondence is the AdS/CFT or gauge/gravity 
duahty in a long wavelength regime. This duality relates a particular 
strongly coupled non-abelian gauge theory in d dimensions to string 
theory, which in a certain regime reduces to classical gravity m d+1 
dimensions. The regimes we will work with are the planar limit 
in the field theory and the long wavelength regime. The latter is 
a small Knudsen number regime, and the CFT can be described 
effectively by a hydrodynamic description. Hence, we have a duality 
between a d-dimensional relativistic fiuid to a classical gravity in 
d+1 dimensions which is just Einstein's equations with a negative 
cosmological constant in an asymptotically Anti de Sitter space- 
time. We can describe this duality holographically and relate to the 
fiuid as "living on the boundary" of the whole d+ 1 space-time, which 
we will refer to as the bulk. Moreover, the small Knudsen number 
suggests an expansion of the fields, from both side of the duality, in a 
derivative expansion. This will allow us to solve Einstein's equations 
order by order. It turns out that parts of Einstein's equations at 
a certain order implement the stress tensor conservation equations 
of the fiuid at a lower order. This statement will be checked in the 
following chapters, but not before we will give the basic set up of 
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the bulk space-time geometry in this chapter. Two reviews on the 
gauge/gravity duahty can be found in |TT1 [T2] . 

3.1 Preliminaries: Schwarzschild black holes in AdSd+i 

We would like to look at the dual description of fluid from the gravity 
perspective. In order to do so, we will have to solve Einstein's 
equations with a negative cosmological constant with a particular 
choice of units {Racis = 1) and we get: 

Eab = Rab — 2^9 ah 2 3 ah = T'^"'^^^'^ = 

^Rab + dgah = Q, R=-d{d+l) (3.1) 

One solution to these equations is just the AdS^+i solution which 
is dual to a vacuum state in the CFT. Another class of solutions is 
described by the "boosted Schwarzschild black branes", 

ds'^ = -2u^dx^dr - r^f (br) u^u^dx^dx" + r'^P^^dx^dx" (3.2) 
with, 



/(r) = l- 


--,, u>^={l,iu'), 


7 = 




X 


^- 


- («f 




[u') = u'Ui 






(3.3) 



This solution is dual to a CFT at a finite temperature T, where 
the velocity m* and the temperature are constants. 



The metrics (3.2) describe the uniform black brane in an asymp- 



totically Anti de Sitter space-time written in the ingoing Eddington- 
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Figure 3.1: Penrose diagram of the uniform black brane and the causal structure 
of the space-times dual to fluid mechanics illustrating the tube structure. The 
dashed line in the second flgure denotes the future event horizon, while the 
shaded tube indicates the region of space-time over which the solution is well 
approximated by a tube of the uniform black brane. (Taken from |11|). 

Finkelstein coordinates^ at temperature T, moving at velocity m*. 

The sets of solutions characterized by d parameters give us d pos- 
sible different solutions. The parameters are just the temperature 
and the d parameters which define the d-velocity, so we have d + 1 



parameters, and because of the normalization (2.2) they are reduced 
to d parameters. 

If we had written the solution in Schwarzschild type coordinate^ 
we could easily have read from the metric the location of the event 
horizon r = 



One can look on the left image of figure 



3.1 



in order 



to understand in a pictorial way the casual structure of space-time. 



sion 
2 



^ These coordinates exhibit a regular solution at the event horizon for all orders of expan- 
The unboosted solution in Schwarzschild type coordinates 

dr2 



ds^ 



-r^f{br)dt^ 



_ , , -I- r Siidx^dx-' 
r^f(br) ^ 



with 5ij is Kronecker's delta. Note that it is not exhibits regularity at the event horizon 
location. 
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3.2 Away from global equilibrium 

Now we would like to consider an out-of-equilibriuni black brane in 
such manner that the dual description will give us hydrodynamics 
of a viscous fluid j2]. In order to do so, we promote the parameters 
to be dependent on the boundary coordinates x'^, but not on the 



radial direction r, and we get from (3.2): 



ds"^ = -2Uf, (x") dx^dr-r^f {b (x") r) m^ (x") u^ (x") dx^dx'^WP^.u (x") dx^dx"" 

(3.4) 
However (3.4), denoted as f^*^*^) [u (x") , h (x"))p| is not a solution 



of Einstein's equations. 

The procedure to satisfy Einstein's equations goes as follow: 

1. Expand the velocity and temperature flelds in a small param- 
eter. 

2. Use the ansatz so that after a while the system relaxes and go 



back to the static solution (3.2). 



3. Calculate the metric order-by-order to satisfy Einstein's equa- 
tions. 

We will do the expansions in the manner of effective fleld theory. 
As we have done in (2.7[), we will take all possible derivatives along 



the boundary coordinates (i.e., x^) to get the expansions, 

oo oo 

«^(x") = 5^«5)(x"), 6(x") = 5^6(o(x") (3.5) 

1=0 1=0 

It is possible to prove that every derivative will come with a 
factor of b, and if the flelds are changing slowly when looking on a 
large length scale L, we have that every term in the expansions is 



Consider the fields as the zeroth order of 3.5 
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Figure 3.2: Cartoon of 'tubewise approximation' of slowly- varying configuration 
by a corresponding piecewise-constant one. (Taken from [12!). 

proportional to(-|) ^ {yl) ■ In other words, in order that these 
expansions be justified we will require the Knudsen number to be 
small 12.61 

This regime of expansions can be visualized as tubes that go 
along the radial direction and their basis is on the boundary (figure 



3.1 right picture). In those tubes we can define a definite velocity 



and definite temperature (figure 3.2). 



Now we will plug the expansions (3.5) into the metric and get 



Qab = Yl'i^o9ab- Then we will require that this metric solves Ein- 
stein's equation and we can determine the metric order-by-order in 
an iterative form. 

The Einstein's equations can be written in the form: 



(3.6) 



where EI is a linear second order differential operator in the r variable 
alone. Note that it does not depend on n, so that it is the same 
operator in any order of the expansion. Moreover, the precise form 
of this operator at the point x'^ depends only on the values of 6(o) 
and mJ'qn at x^ but not on the derivatives of these functions at that 
point. The source term s„ is a regular source term which contains 
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complicated combinations of derivatives of the fields and the fields 
themselves and it is different in any order of expansion. 

It turns out that it is possible to classify Einstein's equations into 
two categories: 

1. the constraint equations E^ = 0. 

2. the dynamical equations E^ = 0, E'^ = 0. 

The d constraint equations are defined as the equations that are 
of first order in r derivatives. It follows that they are connecting 
the derivatives of the temperature field with the derivatives of the 
velocity field. Hence we can use them in order to eliminate the 
temperature derivatives in all order of expansion. The remaining 
^ ^ ' dynamical equations are determined by the next order metric '^ 
Then all together we have -^ which is the number of equations 



that (3.6) possess. 



3.3 The Bulk Metric 

In this section we will present the general form of metrics that can 
describe the bulk geometry of AdS and are also in accord with the 
conformal fluid on the boundary of the AdS space-time. Then we 
will present the second order solution of the metric which satisfy 



(3.1) in a Weyl covariant formulation |13) . 

We start by writing the most general metric that satisfies the 
gauge choice of the previous section: 



ds"^ = -2m^ (x") dxf" {dr + V^ (r, x") dx") + 0^^ (r, x") dx'^dx'', 

(3.7) 



*There is a redundancy among the remaining equations which leaves — Hs — '- independent 
'dynamical' equations, this freedom is eliminating by choosing the following gauge condition: 
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where 0^;^ is transverse to the fluid velocity, i.e., u^&^y = 0. 
We will also present the inverse of the bulk metric. 



U>' [{d^ - V^dr) ® 9, + 9, (d^ - V^dr)] 

or in another form. 



(3. 



(3.9) 
where the symmetric tensor (C5^^)'"^ is uniquely defined by the re- 
lations: u^ {(d-^)"" = and {(3'^y^(3xu = Plf- 

We will now write the solution to the second order. 



V. 



rA 



^ ' d-2 

u 



D,ui - D^a] 



n 



2{d-l 



-u,, 






"M 



2{br) 



r' (l - (br)') - \uj^,u-^ - {brfK, (br) ""-f^ 







liv 



and. 



+2{brf F(br) 



-2 {brf Hi {br) 



-a^^ + F {br) a^axu 



2{brfKi{br)^^^^P, 



d-l 



flU 



yal3 



U L>xO-uu + CTu ^Az. -, —P, 



d- 1 



flV 



+2 (6r)' H2 {br) [u^Dxa^, + u^axu + w>^a] 
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-l\fiU 



(0-) 



.p^- + -w'^^^ 



262 



F{br) 



-a^"" - F (br) a^^a 



A* Xu 



%K, (br) ^^P- 



d- 1 



+—H^ (br) 



u 



a— I 



26^ 



H2 (br) [u^D^a^^" + u^y + u\a^'^] 



with. 



F{br) 
Hi (br) 



00 y,-l 



1 



00 yd-2_i 



br 



yiy'^ -I 



dy 



-dy 



5 .,d~2 



y^~' - 1 



«,(,.)4pw^-/;^)^^ ^,^._„ 



rfy 



iri(6r)= r^ Tdyy'F'iy) 

Jbr S J£ 



i^a (br) 



br 






i-e(e-i)i^'(0-2(rf-i)e 



d-l 



/oo 
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Chapter 4 



Horizon Dynamics 



Given a manifold of space-time one would like to investigate the 
geodesic dynamics equations of that manifold. To do that one has to 
solve the Einstein equations. However, it is possible to consider slices 
of that manifold, for instance hyper-surfaces. Einstein's equations 
projected onto a hyper-surface take a simpler form called the Gauss- 
Codazzi equations. 

In order to understand the geodesic dynamics on hyper-surfaces 
of the space-time geometry (Gauss-Codazzi equations), we need to 
introduce some basic quantities that characterize the hyper-surfaces 
(induced metric, extrinsic curvature, etc.). 



4.1 Geometric Preliminaries 

A hyper-surface is a sub-manifold that can either be space-like, time- 
like, or nulP](we will almost not deal with the latter in this section). 
A particular hyper-surface S is selected either by imposing a restric- 
tion on the coordinates, $ (x*^) = 0, or by providing a parametric 
equation of the form x"' = x°' {y^). 

The normal vector to the surface S is defined by $ „, because the 
value of $ changes only in the direction orthogonal to S. 

^ For the null case refer to [151 116] 
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A unit normal m°- is defined by: 

$ 

rua = , '" (4.1) 

where we demand that m" point in the direction of increasing $: 
m"$ a > 0, and we get its norm to be: 

, -1 if S is space — like 
rrfrria = { (4.2) 

+ 1 if S is time — like 

In the case of a null surface we can only define a class of null 
normal co-vectors, in the following way: 

L = f^,a (4.3) 

where / is an arbitrary function. The norm of the class of vectors 
is: 



(4.4) 



Now we would like to define a metric that would be intrinsic to 
the hyper-surface. The induced metric h^,, to a hyper-surface S is 
obtained by restricting the line element to displacements confined 
to the hyper-surface. Using the parametric equations x"' = x"" (y^), 
we find that the vectors, 

are tangent to curves contained in S. (This means that e^m^ = 
Cp|. Therefore, we now write the restriction of the line element for 
displacements within S: 

■^e"ma = oc TTu TT-a = TT-u = 0. It is zero because the function is constant in the 
directions of the surface. 
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ds^ = Qabdx^'dx 






where 



V = dabe^.et (4.7) 

Note that the induced metric is a scalar under bulk space-time 
coordinate transformations, a;" — )■ a;" . However, it transforms as 
a tensor under the hyper-surface coordinate transformations, y^ — ?■ 

With the aid of the induced metric we can write the decomposi- 
tion of the bulk metric: 

gab ^ j^ab ^ ^a^b ^4 g^, 

with: 

^ab ^ /^"/3e»eg (4.9) 



where the plus sign in (4.8) represents the time-like hyper-surfaces, 
and the minus represents the space-like hyper-surfaces. 

This decomposition can be defined to every hyper-surface that is 
not null. 

We now would like to understand how tangent tensor fields are 
differentiated. First we define such fields: 

A tangent tensor field A"-^'" is a tensor field that is defined only 
on S and is purely tangent to the hyper-surface. Such fields admit 
the following decomposition: 
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^ab... _ ^/.^...gagfe . . . (4.10) 



Equation (4.10) implies that A°'^"'ma = A"'^"'7hi, = ■ ■ ■ = which 
confirm that A""^'" is tangent to the hyper-surface. The tangent 
tensor A^*^ indices are lowered and raised with the induced metric 
on the hyper-surface h^i, and h'^". 

Now we can define an intrinsic covariant derivative. We will 



use a tangent vector field for simplicity: A"- = A^e'f^., A^-fha 
0, A^ = AaC't- The intrinsic covariant derivative of A^ is defin 
by the projection of Aa-^i, onto the hyper-surface: 



A,\, ^ A,,,elel (4.11) 

It is possible to prove that equation (4.11) can take the form of, 



A,\. ^ A^,. - Tl^A^ (4.12) 

where 

^% = ^"^ (V,- + ^-,/. - ^M--.) (4-13) 

We can now introduce the extrinsic curvature of a hyper- surface, 
which is just the projection of the covariant derivative of the normal 
to the hyper-surface onto the hyper-surface S: 



K^, = rha-be^^el (4.14) 

This quantity is a symmetric tensor and it tells us how the hyper- 
surface S is embedded in the bulk. We can see this by taking its 
trace: K = h^^K^j,i, = m"^ (we used the decomposition equation 
(4.8)) which is the expansioEp] of a congruence of geodesies that 



intersect the hyper-surface 5* orthogonally. When i^ > (the con- 



^This term will be explained in chapter 
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gruence is diverging), the hyper-surface is convex and if i^ < (the 
congruence is converging), the hyper-surface is concave. 

4.2 Gauss-Codazzi equations 

We are ready to introduce the Einstein equations on a hyper-surface 
S. We start by defining a purely intrinsic curvature tensor by the 
relation: 

Al. - A\u, = RL,A' (4.15) 

where the explicit expression for the purely intrinsic curvature tensor 
is: 

Rxu^l = ^x^i,u ~ ^\u,^l + ^au^x^J. ~ ^aiAxu (4.16) 

We now want to establish a connection between the intrinsic cur- 
vature tensor to the Rieniann curvature tensor defined in the bulk 
geometry of space-time. It can be proven [14j that the relation is: 

(4.17) 
where again the plus sign represents the time-like hyper-surfaces and 
the minus represents the space-like hyper-surfaces. 

Projecting along m^ and using the symmetry of Riemann curva- 
ture tensoij^ gives: 

Radcbm'^ele^f^e'^ = {K^p^^ - K^^i,s) (4.18) 

We will now contract the indices by using the metric h"''^ with 



equation (4.18). Then, by using equation (4.8) and the symmetry 



of Riemann curvature tensor, we get: 



23 



Ram^e'p = {K1\^ - K^p) 



(4.19) 



If we look at the bulk metric Qab contracted with ■nf'e\ we get, 



Qabrrfe^p 



rhbC^p 



(4.20) 



This means that we can replace the left hand side of (4.19) with 



{Rab + Cgab) rh^-e^a, where C is an arbitrary coefficient, because every 
term that is proportional to the metric will vanish. For instance, 
we can replace the left hand side with the Einstein tensor Gat = 
Rab — ^Rgab to get the Gauss- Codazzi equations, 



GabTh^ = i^||„ - K^p (4.21) 

Because we are solving the Einstein equations with a cosmological 



constant we would like to replace the left hand side of (4.19) with 
Eab, which we introduced in equation (3.1), and we get. 



Eabm'^e' 






K 



,/3 



(4.22) 



As mentioned in the equations (4.21 ) and (4.22) we are eventually 



reducing these equations, because of (4.20), to (4.19) and these are 



the equations we have to solve in order to find the geodesic dynamic 
equations. 



4.3 The Constraint equations 

In the previous chapter we stated and showed that Einstein's equa- 
tions can be divided into two categories, the dynamical equations 
and the constraint equations. In this section we will deal with the 
latter. 



Claim : The constraint equations E'^ = Q with the metric (3.7), 
which solves the dynamic Einstein's equation, can be written in the 
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form: 



Eab^el = 



(4.23) 



where ^'^ is normal to any hyper-surface 5* (time-hke, space-hke and 
null) that is defined by $ = r — r (x") = 0. Thus ^^ is defined by 
^^ = f g^°'^a, where / is an arbitrary normalization function. 

Proof : Let us set / = 1 for convenience reasons and look at the 



LHSof (4.23): 




Eab^ e^ = Rab^ e° 




= e^^Rabg'^d,^ 




= elK-elR:d,r{x-) 




= R; + R;d,r{xn-R';d,r{x'^) 


-K5^r(x")9,r(x° 



,^ ^„^, ,^ , , ^.^^,. ,^ , {g''''R„r + 9''''Rrr)d,r{x'')d,r{x'') 



R^ - dd^,r (x") + dd^^d^r (x" 
R^ + dg'^^g^rd^r (x") d,r (x° 
""u^df^r (x") d,yr (x" 



r:. 



where we got the fifth line by using Einstein's equations (3.1) and 
we got the sixth line by using the relation: g^^grv = —^^'^u^ = 0. 
Using the relation: E^ = i?^ completes the proof. 



For the non-null case we can take ^° = m'j^in (4.23) and we get 



from (4.22) the Gauss-Codazzi equations with a vanishing matter 



stress tensor: 



^fi\u ~ -^,M 



(4.24) 



'Note that the normalization function / can be different from one. 
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If we look at the constraint equations on the boundary with the 
normal vector to the boundary surface ^°' = g"-"^, then the constraint 
equations, according to jlTj, take the simple form of the equations 



of motion for the boundary fluid (2.5). 



The constraint equations just connect the derivatives of the fluid 
d-velocity with the derivatives of the temperature (or the inverse 
temperature b), and of course the rf- velocity and the temperature 



fields depend, by our definition (3.5), on the coordinates perpendic- 



ular to the r coordinate, which are the boundary coordinates x". 
For this reason we can claim that the constraint equations are the 
same for every hyper-surface ^ = r — r (x") of space-time. In the 
following chapters we check this claim explicitly. 
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Chapter 5 



The Second Order Constraint Equations 



5.1 Hyper- Surfaces At Constant Radial Location 

In this section we will check the constraint equations to the first 
order in derivative expansion for a constant radial coordinate hyper- 
surface, i.e., ^ = r — R, where R is constant. 

For reasons of convenience, we will write here only the zeroth and 



first order parts of the metric (3.7) 



ds'^ = -2Uf,dx^dr-r^f{br)u^u^dx^dx'' + r'^Pf,udx^dx'' (5.1^ 

2 
+2r%F (br) a^^dx^dx" + - — -ru^u^dxu^dx''dx'' - ru^dx (u^u^) dx'^dx'' 

The components of the inverse metric to the first order are: 



a — I r J.Z J./ 

(5.2) 



/^ = r^f(hr)-- -rdxu^, g'"" = u^ — m^^aM^, ^^" = -P'^" ^^a''" 



Our basis vectors (4.5) are: 



dx^ ~ ' ^~ dx^" 



< = 7^ = 0, el = ^—=5l (5.3) 



27 



The normal vector to this surface is: 



m 



g^^db^ 



y/\g->'da<^db<^\ v1? 



(5.4) 



The components of the one-form dual to the normal vector are: 



rrir 



•/\ 



m. 







9"^ 



(5.5) 



In order that the last expression will be defined, we demand g^"^ ^ 
0. The surface $ is time-like if g" > 0, or space-like if g^"^ < 0. 



Our goal is to write the constraint equations by using (4.24), so 



we will need to find the extrinsic curvature of the hyper- surface. 
We will calculate the mixed upper and lower indices of the extrinsic 
curvature in the following way: 



Ki^ 



h^^ma;bel 

{g^''Tm^m'')ma;bel 
9^''rha-bel T m^m''ma;bel 

'. u b 



(5.6) 



where we used (5.3) to get the fourth line and m°'rha\b = to get 



the sixth line. After applying the covariant derivative, we get: 



ir^ = V^m'^ = d^m^" + V^Jh'' = d^m" + V^^m" + V^^rh" (5.7) 

In order to compute the constraint equations up to the second 
order in derivative expansion, we have to find the constraint equa- 
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tions that come from the zeroth order metric, and then apply them 
to the calculation of the constraint equations that comes from the 
first order metric. We show here the explicit form of the extrinsic 
curvature to the zeroth order, 



_^/^(o) 



and its trace is. 



1 



1 d 



V\TWY\\ 2(6ry 



juf^Uu + f (for) T]^ 



(5.8) 



Kio) 



1 d 



VWWTlVibr) 



+ df{br) 



(5.9) 



Inserting (5.8) and (5.9) to the Gauss-Codazzi equations (4.24) 
we get: 






d 



{2duyU^d^ In 5 - u^d^u^ 



-u„d^u^ + (d + 1) d^ In h 
-dPl^df,\nb) = (5.10) 



Projecting (5.10) along u'^ we get: 



2 {bRf vwJm 

^ d^u" - {d - I) u^d^r^b = 



id^u^'-id-l)u''d^\nb) = 



(5.11) 



Projecting (5.10) along P^ we get: 



d 



2 (bR)" .ATJm 

> uf'd^u^ - P^d^lnb = 



uf^d^u^ - P^d.lnb) = 



(5.12) 
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Now we compute the extrinsic curvature to the first order from 



the first order bulk metric (5.1) 
_^m(o+i) ^ ... 



, , — r :iU^Uy + rf (br) ri^, + o"'^ 



-r^f (br) drF (br) ba^i 
2 / (br) (br) 



(5.13) 



j^{0+l) ^ ^(0) 



(5.14) 



We define the quantities with the superscript: *^°\ '^^\ ^^^ to be 
quantities that have no derivative, only one derivative or two deriva- 



tive terms respectively. Now we insert (5.13) and (5.14) to Gauss- 



Codazzi equations (4.24) and get the constraint equations to the 
second order: 



(0+1) _ (0+1) 



_f'^(l);^M(l) _ f o-(2) 7>'/^(0) 

^y (7 flV (T 

2(bR)''^\f(bR)\ d ^ ^ "" ^ ' ^ "' 

(5.15) 



where the bar over the Christoffel symbols represents the Christoffel 
symbols which have been calculated with respect to the metric h^y 



of the hyper-surface $. The fifth line is the LHS of equation (5.10). 



Projecting (5.15) along u'^ we get, 
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, '^, Id.u^' -(d-1) u^d. In 6 - ^(T„fl(T"^ 

2{hRY,/\nm\ ^ ^ ^ d "-^ 

2h 
^ d^u^" -(d-l) u^d^ Info = —aa(s(T''^ 



= 0^ 
(5.16) 



where we used the relation, 



-d.ay' 









d 



O'a/SO- 



aP 



Projecting (5.15) along P^ 



2 {hRf VlTIM 

X (-u^d^u, + P^d, In h + ^P; (a^< -{d-l) A^a^^) ) = 

^ u^d^u^ - P^d^lnb = ^P; {d,a^ -{d-l) A^<) (5.17) 

We can see that our results are in agreement with the equations 



of motion (2.5) of a conformal fluid to the second order in derivative 



expansion, that is, up to first order terms of the stress tensor (2.10) 



5.2 Non-Constant Radial Location 

We can generalize our discussion to non-constant hyper-surfaces, i.e.: 
R = R (x). In order to do so, we need to use the property in which 



31 



we can explicitly calculate the r derivatives of tensors which have 
been calculated in the bulk. For instance, the constraint equations, 
E"^, = i?y, = can have the followiner form: 



R"^ = y^ Qi (r, b) hi^ {du (x) , dh (x) , ddu (x) , ddh (x) , higher terms) 

(5.18) 

Note that in this expression there are no derivatives of r because 

we can calculate them explicitly. However, derivatives of x^ can 

not be calculated explicitly, because they are operating on arbitrary 

functions u^ (x) , 6 (x). 

For the second order constraint equations, we have no higher 



terms in the brackets of (5.18) and the answer for constant R was 



calculated above by Gauss-Codazzi equations (5.15). From (4.19) 



and the claim in chapter 4 (4.23) we have: 



jLpr(l+2) 



KT'^-Kl'^'^=9i{R)h. = (5.19) 



where: 



91 (R) 



d 



2 (bR)' VWWWl 



(5.20) 



hu 



2dui,u^dfj_ In 6 — u^d^u^^ 



Uyd^ + {d + I) duhih - dPj^df, In b 
2h 

-(9X-(rf-i)^X) 



(5.21) 



and 
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A^ = ^/\g''^^a {r - r (x")) db{r-r {x»))\ (5.22) 



Therefore, the sum in (5.18) has only one term: gi (R) = Vz > 2. 
However, the latter is supposed to be true for \/R = const. Thus: 
Qi (r) = Vi > 2. Therefore, we obtain the constraint equations for 
the second order in a derivative expansioiiQ which are the same for 
any non-null hyper- surf ac^ 



^Because gi (R) ^ 0, we can divide the equations: gi (i?) hi^ = by gi (/?). 
■^However, the constraint equations for / {bR) = 0, which is in the case of a null hyper- 
surface, hold true and are still the same as in the non-null case i6j. 
The case in which ij = is not define. 
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Chapter 6 



Event vs Apparent Horizon 



In this chapter we review the different notions of black hole horizons. 
We will focus our discussion on two distinct horizons, namely, the 
apparent and the event horizons. The last section in this chapter 
is devoted to calculating the constraint equations from the second 
order bulk metric projected on the different horizons. 

6.1 Event Horizon 

One approach to define the unique hyper-surface that defines the 
boundary of a black hole is the causal approach. In this approach one 
defines the unique hyper-surface, which is called the event horizon, 
as the boundary of the casual past of future null infinity [18]. This 
mathematical definition can be put in a more physical language by 
the following statement: the event horizon is the boundary of the 
region in space-time from which nothing can ever escape. 

We should note here that the future null infinity is a time-like 
hyper-surface. That is because we deal with an asymptotically AdS 
space-time, for which the future directed null geodesies (which is 
the future null infinity in the definition) end on a time-like hyper- 
surface. 
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One should also note that by definition the event horizon is a null 
hyper-surface which is generated by null geodesies. 

We would like to describe a method for finding the location of 
the event horizon. This method requires the knowledge of the late 
time generators of the event horizon, which means the location of 
the event horizon in late time. The method goes as follow: evolve 
the geodesic equation backwards in time and use the knowledge of 
the late time generators as the future boundary condition. 

For our metric, which describes a dynamical black hole which is 
dual to viscous fiuid fiow, we can assume that after a while the fiuid 
relaxes, which means that we have a stationary black hole that is 



described by the metric (3.2). The location of the event horizon for 
the metric is well known and is r'j!^'^ = ^. Since we are working in 
a small Knudsen number i.e., derivative expansion, we can find the 
location of the event horizon order by order |17) : 

oo 

rEH{xn = Y.'-EHi^") (6.1) 

1=0 

We require that the normal vector £" to the event horizon hyper- 
surface Seh , defined by: 

Seh (r, x) = r- teh (x) , t = g'^'d, Seh (6.2) 

will be a null vector i"-ia = and we get the equation, 

g"- - 2d,rEHg'^ + d^rEHd^rEHg"" = (6.3) 

which can be solved algebraically by using the initial condition for 
the location of the stationary event horizon: r^^ = ^ . The solution 
up to the second order. 
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Figure 6.1: The event horizon r — rEnix'^) is sketched as a function of the time 
t and one of the spatial coordinates x (the other d — 2 spatial coordinates are 
suppressed). (Taken from jll)). 



TehKX 
with 



h (x) 



+ h {x) {hiaa^pa''^ + /i2C^a/3C^"^ + h^n) (6.4) 



hi 
h2 



2{<f + d- 


4) 


d^ {d - 1) {d 
d + 2 


-2) 


2d{d-2) 
1 





^2(1) 
d{d-l' 



d{d- l){d-2) 



(6.5) 



Note that the solution is dependent on the boundary coordinates 
alone, which means that for a constant x'^, i.e., a tube along the 
radial direction, we have a fixed solution and that in late time we 



get our starting condition that is the static solution (figure 6.1). 
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6.2 Apparent Horizon 

We would now like to describe a black hole by relating to its strong 
gravitational fields instead of its causal structure. The latter defini- 
tion led us to the introduction of the event horizon hyper-surface in 
the previous section. In this section we show the major drawback in 
that definition, which leads us to a new hyper-surface, the apparent 
horizon |19| . 

In order to find the location of the event horizon, one needs to 
know the entire future evolution of space-time, which means that 
one needs to wait infinite time in order to find the location of the 
event horizon. To illustrate the problem, we will present the ex- 
ample given in [20j. Imagine a Schwarzschild black hole which has 



been irradiated by an infalling shell of null dust (figure 6.2). This 
is represented by the gray rectangle. The inward and outward null 
directions are tangent to the light gray dashed lines. Note that the 
inward null direction is horizontal. In order to find the event hori- 
zon, one needs to track the evolution of radial null geodesies. If they 
fall into the singularity, they are inside the black hole and if they 
are heading outward, even after the null dust passes, then they are 
outside the black hole. If one wants to provide the location of the 
event horizon without having any prior knowledge of the arrival of 
the null dust, one would state that it is located at the hard dashed 
line, which is clearly a false statement. One can realize this by ex- 
amining the radial null geodesic B, which seems to be escaping from 
the black hole, but ends at the singularity. By a close examination 
one sees that the location of the event horizon is the hard line. One 
can see how it curves a bit in anticipation for the null dust to arrive, 
as if it "knows" that it is coming. 

We will now explain how a new perspective on a dynamical black 
hole avoids the causal problem. This perspective is, as mentioned 
before, the strong gravitational field approach. In order to describe 
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Figure 6.2: A schematic demonstrating the non-local nature of event horizon 
evolution for a spherically symmetric space-time with the angular dimensions 
suppressed. Horizontal location measures the radius of the associated spherical 
shell while time is (roughly) vertical. The shaded gray region represents infalling 
null dust. (Taken from pO]). 



this approach we will give an example found in [T^ I2U] . Imagine a 
spherical shell, in a four dimensional world, which emits light and 
then stops. The light front null geodesies will be divided into two 
families that move in two null directions: the outgoing null direction 
and the ingoing null direction. We will denote the tangent to the 
light ray null vectors, which are perpendicular to the spherical shell 
surface, by f^ and n"' respectively. The outgoing light front will 
increase in cross-sectional area and diverge, while the ingoing light 
front will decrease in the cross-sectional area and converge. In a 
more mathematical way, if we define the expansion parameter 6 |14| 
as the fractional rate of change (per unit affine-parameter distance) 
of the congruence's (a family of geodesies that does not intersect 
each other) cross-sectional area: 



e 



1 d 

JaUx 



5A 



(6.6) 
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trapped^surfaces 

apparent horizon 




Figure 6.3: An "instant" Sf along with some of its trapped surfaces (small black 
circles), the associated trapped region (dark gray) and the apparent horizon 
(thick dashed line). (Taken from j20]). 

where A is an affine-parameter, and 6A is the infinitesimal cross- 
sectional area measured in the purely transverse direction^ then 
we can describe the convergence and divergence of the light front by 
the following: 

^(,) > 0, ^(,) < (6.7) 

where the sub-script implies to which light front the expansion is re- 
ferring to. A positive expansion means that the cross-sectional area 
will increase while negative expansion means that it will decrease. 

Now we are ready to introduce the notion of a trapped surface: 
A trapped surface in d + 1 space-time dimensions is a space-like 
hyper-surface oi d — 1 dimensions, which, by definition, has two null 
normals: the outgoing null normal £° and the ingoing null normal 
n", where the expansion in both directions is negative. 

^(,) < 0, 0(„) < (6.8) 

The hyper-surface that surrounds all of the trapped surfaces is 
called the apparent /ionzorip] (figure 6.3) and defined mathematically 
by: 



^The geodesies are crossing the cross-sectional area SA orthogonally. 

■^The apparent horizon is a special case of the so called marginally trapped surface reviewed 
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^(,) = 0, ^(„,) < (6.9) 

Note to the condition 6'(^) = 0. It does not mean that if matter 
comes into the black-hole the apparent horizon location will not 
change (increased If we let the apparent horizon evolve in time, 
we get a co-dimension one tube (which means the dimensions of 
the tube is d). The tube is built in the following manner: given a 
foliation of space-time with d— 1 hyper-surfaces Sx which have two 
null normal vectors: i"" and n"' as defined above, we connect them 
with an evolution vector field z/" that maps every point in ^a to Sa+i 



to create the c?- dimensional tube A (figure 6.4). The evolution vector 
is perpendicular to each of the hyper-surface ^a and tangent to the 
d-dimensional tube A. If the Sx are the apparent horizon hyper- 
surfaces, then we say thatA is a time evolved apparent horizor^ 

We will now explain about the geometric quantities that we de- 
fined. First we state that we cross normalized the null normal vec- 
tors, 

tua = -1 (6.10) 

Moreover, we present their normalization: 

r ^ /r n" -^ ^n"" (6.11) 

where / is an arbitrary function. 

Then we define the evolution vector to be normal to the ^a hyper- 
surfaces, so it will be built from the two null vector^ f^ and n": 

•^The expansion of the apparent horizon is related the expansion in the direction of the 
evolution vector i/" ( |6.12[ l, which is: 

^{i/) = ^(i) - CB^^n) = -C9(„) > 

*Froni now on we will refer to the time evolved apparent horizon as the apparent horizon. 
^Because we are using null vectors, we can always choose a normalization that will give us 
this form of evolution vector. 
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Figure 6.4: A schematic of an d-tube A with compact foHation hyper-surfaces 
Sx along with the outward and inward pointing null normals to those hyper- 
surfaces, ly'^ is the future-pointing tangent to A that is simultaneously normal 
to the Sx. (Taken from [20,1. 



-Cn"" 



(6.12) 



for some function C Note that if C > 0, then we get space-hke 
hyper-surface, if C = 0, we get a null hyper-surface and if C < 0, 
we get a time-like hyper-surface. 

Now we can present m", the normal to the d-tube hyper-surface 
A vector, which is perpendicular to the evolution vector and we get: 



m 



r + Cn" 



(6.13) 



Because m"' is normal to the hyper-surface A, which can be pre- 
sented by the hyper-surface equation (p (r, x) = 0, we will choose to 
define it as: 



m» = g^^db(f 



(6.14) 



which is similar to what we did in (4.1) but without having to nor- 
malize itPl 



^Note that m" can be space-like, time-like and even null. 
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We pause here and return to the event horizon description. The 
event horizon can also be described by the rf-tube description and 
of course its evolution vector will be with C = 0, so we actually 
see that both the normal vector to the hyper-surface A (the event 
horizon) and the evolution vector are the same: 

m'^ = i/'^ = t (6.15) 



where i"- is defined in (4.3)''' 



We now specify some more geometric quantities and we introduce 
the metric for the hyper-surface ^a, as well: 

Qab = 9ab + Lnh + hria (6.16) 

Note that this metric has two zero eigenvalues. 
We can associate the expansion parameters to the geometric 
quantities using the following definition: 

^(£) = g'^'VJb, ^(n) = g'^'Van, (6.17) 



Note that q""^ is defined as in (6.16), but with upper indices and 
it is not the inverse of Qa^ 

Now we are ready to return to the causal problem (the global 
definition of the event horizon) presented earlier and see whether 



we can provide a solution to it. One can see (figure 6.5) that the 
apparent horizon changes only when the null dust comes in, as op- 
posed to the event horizon that "anticipates" the arrival of the null 
dust, and curves even before it arrives. Therefore, we have found a 
solution to the problem. Note that the apparent horizon is located 
inside the black hole, which means that it lies within the event hori- 
zon. Moreover, the apparent and the event horizons are the same 



^In the rest of the section we wiU set the normahzation of / in 1 4.3 i to one. 

*One can think of q°-^ as a projection onto the S\ hyper-surfaces, which in | |6.17| restricts 
Va^6 ^nd Va"-!) to the hyper-surfaces. The expansion 6 is the trace after the restriction. 
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Figure 6.5: A simulation similar to that of Figure 4 though this time two distinct 
shells fall into the black hole. Both the apparent and event horizons are plotted. 
(Taken from j^P]). 



in the stationary case; one can see this by looking at their location 



after the two null dusts pass (figure 6.5). 

We want to find the location of the apparent horizon in the same 
manner that we did for the event horizon, by using a derivative 
expansion: 



rAuV^ 



E 



,(0 

AH 



\X 



(6.18) 



From the requirement 0(^) = we can find the location order by 
order. 

However, in order to do so we need to foliate space-time with 
d — \ hyper-surfaces S\. We define those surfaces by specifying 
the two null vectors to the surfaces. The symmetry of space-time 
singles out the ^a hyper-surfaces of constant r and constant "time" 
propagation. The latter is an arbitrary combination of x^ . The 
outward null normal 1°" and the inward null normal n°- are chosen to 
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be: 



dr J \ ^x^^ 



n 



dr 



(6.19) 



(6.20) 



We can determined them by the relation (6.10) and the fact that 
they are both null vectors and we get 



AT^ = M^ + r^, r^u. = 



(6.21) 



B 



-M^V, 



(6.22) 



to the second order, where T^ is composed from one or two deriva- 
tives of the d-velocity or the temperature, so we have freedom in 
determining this quantity. We therefore require that this descrip- 
tion be suited to describe the event horizon, because we know that 
the apparent horizon coincides with the event horizon in the station- 
ary case. In order to do so, we require that the null outward vector 
in the x^ direction resembles the null outward vector in the x^ di- 
rection that defines the event horizon. This requirement determines 
T'^, and we get: 



r^ 



1 - 1) bu'd.u^ - ^^P- (9a -{d-3) A,) u^i + 



+ '^r(^^) + Z^ 



2b' 



^ lP^''{dx-{d-l)Ax)a^p + 



+ I 1 - ^"j ^F (6r) a'^-'u^dxH, 



(6.23) 



Now we are ready to find the expansion in the outward null di- 
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rection: 



\i) 






+ 



2r (6r) 



r" {\ - {hrf\ 



2(brf 1 



a/3_ 



d-2 



d-2 



K2 ibr) , 2 a/3 



r^ f 1 — (6r) 
2 (6r)'' 



+ 



7^ 



2r (d - 2) 



-fe^S^iTi (6r) (T„;3a"'5 + -^oj^puj"^ 16.24) 



Inserting (6.18) to (6.24) and using d{^t) = 0, we find order by 
order the location of the apparent horizon: 



rAH[x 



b{x) 



+ b (x) (hicr^^cr''^ + /i2Wa/3w"^ + hsV) (6.25) 



hi = hi 



d? {d - 1) 



(6.26) 



where hi, /i2, h^ are defined in (6.5) 



Note that up to the first order in derivative expansion, the event 
and apparent horizons coincide. 

6.3 Event and Apparent Horizons Constraint Equations 

We are now well equipped to check if the constraint equations for the 
event and apparent horizons are the same up to the third order in 
a derivative expansion. The work that was done in JU\ showed that 
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the event horizon, to the first order in derivative expansion of the 
metric, exhibits the same constraint equations to the second order in 



derivative expansio#l as in the case of non-null hyper-surface (5.16) 



(5.17). As we mentioned above, the event and apparent horizons 
differ only in the second order in derivative expansion. Therefore, 
the constraint equations up to the second order are the same in all 
of the hyper-surfaces mentioned above. Recall our statement about 
the ability to get the same constraint equations from different hyper- 
surface^^ we see that it turns out to be right up to first order in 
the metric. 

In order to check this statement for the event and apparent hori- 
zons to the second order in derivative expansion of the metric, let 
us look directly at the constraint equations i.e.: EJ^ and calculate 
their difference while they are evaluated on the event and apparent 
horizons hyper-surfaces: 



EH -AH = E'l -En = Rl\ - i?M (6.27) 

f^'r=rEH >^'r=rAH f^'r=rEH f^'r=i^AH ' 

Bear in mind that the difference in W,, between the event and 
apparent horizons to the third order will only come from the eval- 
uation of i?y, in the bulk to the first order and from the location of 
the horizons to the second order. Recall that the difference between 
the horizons location is only in the second order. 

Now we present the calculation of the tensor W to the first order 
in the bulk: 



^T = \^'^-d ("^M In h + du^dp In h - u^dau" - u'^d.u^) (6.28) 
2 (6r) 



^The constraint equations always have one derivative more than the derivative expansion 

of the metric 

10/ 



-^Actually for hyper-surfaces of the form (f> = r — r (x"). 
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From Einstein's equations (3.1) we know that R^ 



jKi) 



for any 
r. Therefore, the expression in the brackets in (|6.28|) is zero. This 



expression is just the constraint equations we found in (5.11) and 



(5.12) 



Now we are ready to calculate (6.27) to the third order: 



EH -AH = i?;W| -K/^^l 

^ 'r=rEH f^ 'r=rAH 

d 1 



d 



2 {brEHY 



-2^^^(^ 



{d^lnb 



H 



+du^dp\nb — u^dau" — u^dyUf^) 







(6.29) 



where we get the fourth line by using the constraint equation to the 
first order. 

There is another method to get this result. This method is more 
cumbersome than what we presented here, however, it poses chal- 
lenging computational techniques; therefore, we presented it in the 
appendix. 
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Chapter 7 



Hydrodynamic Quantities 



7.1 Entropy Current 

In [17j the authors introduce a boundary entropy current obtained 
from the event horizon by using the area of the black hole event 
horizon. By using the connection between the area of a hyper- 
surface to the square root of the metric determinant of that hyper- 
surface, and by using a pull-back along ingoing null geodesies to that 
area, we get the entropy current on the boundarj^} In order to do so, 
we split the boundary coordinate x^ to (z/, x*) and we continue to use 
the same coordinates on the hyper-surface $ = r—r (x). We will get, 
as in the d-tnhe description, a space-like hyper-surfaces S,^ which is 
a constant u slice. This, as well, will divide the components of m^ 
into (r7i'',r7i*)r| After we split space-time, we define the metric on 
the space-like hyper-surface S,^ by the restriction of the line element 
to be only on that hyper-surface: 

ds'^ = gabdx"'dx^ = qijdx^dx^ (7-1) 

^Of course not every hyper-surface is adequate. We require that the hyper-surface admits 
the area increase theorem, that gives us a natural candidate: the event horizon. However, it 
is not the only hyper-surface with that property. 

■^Here u does not represent d different coordinates, but a specific coordinate. These nota- 
tions were also used in |17| . 



Then by using the procedure done in |17|we get: 



J'' 



qm^ 



771" 



(7.2) 



where q is the determinant of qij. It can be shown that this form 
of entropy current is indeed covariant [T7]. Also the authors in [17] 
found the entropy current for the event horizon hyper-surface, and 
[T3] reproduced it to arbitrary li-dimensions: 
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with, 
Ai -- 



d-l 



1^ 



+ h''[BiDx(y^'^ + B2Dxu^'''] (7.3) 






-2A, 



Bo 



1 



d{d-2) 



d-2' 
(7.4) 



In order to evaluate the the entropy current to the second order 
in derivative expansion for the apparent horizon, we calculate the 
difference between the two entropy currents: 



TfJ- _ TM 



qm^' 



4 m" 



qrri' 



4 m" 



r=rEH 



(7.5) 



r=rAH 



Because m^ for the event horizon and apparent horizon is the 



same (C.5), we need to evaluate the difference in Wq on the horizons 



The difference to the order we are seeking will come from inputting 
the location of the horizon to the second order into the zero order 
of the metric determinant, 



^0) 



„d-l 



(7.6) 
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and we get the difference between the entropy currents to be: 



7M _ 7M 



1 m^ 
Am" 



r=rEH 



r=rAH 



^d-ll 



„d-l\ 



Amy \ 'r=rEH 

(d — 1) m'^ 



I r=rAH 



.(2)1 



m 



lj-(d~2)/(2) _^ 

M l)-(d-3) 



-CTaiSCr 



aP 



(7.7) 



my d? 
Therefore, we see that the entropy current for the apparent hori- 



zon IS, 



40 J^^ 



+ 6^ [B^D^a^^^ + ^sDaw^^] 



+ 



(7. 



with. 



Ai 



Ai- 



d^ 



(7.9) 



and the other coefficients are the same as the coefficients of the event 



horizon entropy current (7.4). 



We now have to check if this entropy current be suitable to de- 
scribe the notion of entropy, which means we need to answer the 
question: does the entropy current never decreases? 

In order to answer this question [T7] calculated the divergence 
of the entropy current and required the non-negativeness of the ex- 
pression: 



D^J^ > 



(7.10) 
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They did it with arbitrary coefficients Ai, A2, A3, Bi, B2 and 



found that in order that (7.10) apphes, the following connection 
should be valid: 

^1 = 2^3 (7.11) 

This connection, of course, is holds true for the event horizon 
and also for the apparent horizon. Therefore, we can say that it is 
possible to ascribe a different entropy current to the boundary by 
the apparent horizon hyper- surface. 

7.2 Temperature 

Besides the entropy current, we can relate another thermodynamic 
quantity of the fluid to gravity, that is the temperature. 

The temperature of the black hole can be related to the surface 
gravity of the event horizon by Teh = ^- We use a generalized 
deflnition for the surface gravity shown in [20j : 



K^ = -ni^u^Vai^ = -UktVat + Cni,n''VJ^ = ki + Ck^ (7.12) 

to calculate the surface gravity for different hyper-surfaces. 

The result for a hyper-surface located at r = ^ -|- r^'^^ (r*^) has 
second order derivatives) is: 
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d _ dxu^ 
2b ~ d-1 
d 



{d - 3) r(2) + dr 



+ {brfK2{br'''"^'' 



2 (for 

al3 



-iOaP^ 



d- 1 



(7.13) 



and by inserting the location of the event horizon or the apparent 
horizon, we can get their associated surface gravity respectively. 
These quantities may relate to different notion of temperature fields 
by dividing them by 2tt. 

In order to set the surface gravity to be equal, up to 27r, to the 
temperature field of the conformal fluid we set the normalization of 
(6.11) in such a way that we get: Tjiuid = ^- The surface gravity 



transforms in the following way: 



Ko 



/|.,+ ^ln/ 



(7.14) 



where / is the normalization in (6.11), /t^ is the new surface gravity 



and A is an arbitrary parameteij^ 

Under this change of the normalization /, we need to check how 
the other quantities that we defined transform and if this transfor- 
mation changes our results regarding the location of the apparent 
horizon and the entropy current. 

The location of the apparent horizon depends on how the expan- 
sion parameter transforms: 



^U) - f ^{i) 



(7.15) 



'Note that this is the transformation of the change of the parameter A, because if A — > A* , 



then £*" : 



—E'^ = f L 
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and we can see that it does not change the location of the apparent 
horizon. We can see this by trying to find the "new" location of the 
apparent horizon by setting ^|^n = 0. It means that if / 7^ 0, which 
is the trivial case, then 9(^q = which is exactly the "old" location 
of the apparent horizon. 



By the definition of the entropy current (7.2) one can see that it 



is invariant to the normalization transformation. 
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Chapter 8 



Conclusions 



In this thesis we presented a method and different techniques to 
evaluate properties of an out-of-equihbrium black hole that has a 
dual description as a conformal out-of-equilibrium relativistic fluid, 
which "lives" on the boundary of AdS space-time. We set the stage 
and presented the derivative expansion for both the relativistic fluid 
and for the black-hole. We presented a way to compute the Ein- 
stein constraint equations for different hyper-surfaces. The con- 
straint equations turned out to be just the conservation equations 
of the fluid. We found, as expected, that from any hyper-surface 
which is located at r = r (x), we get the same conservation equa- 
tions. This implies that maybe one can describe a hydrodynamic 
duality, not just from the quantities that are related to a certain 
hyper- surface, namely the event horizon, but also to quantities re- 
lated to other hyper-surfaces. We suggested the apparent horizon 
as a possible candidate, and calculated from the apparent horizon 
area the entropy current and from its surface gravity the tempera- 
ture for a conformal fluid. Our work poses a few future questions. 
The flrst is: which hyper-surface describes correctly the conformal 
fluid properties? The second is: how can we use these results for 
the description of the QGP (Quark Gluon Plasma)? 
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Our solution for the difference between the location and the en- 
tropy current of the event and apparent horizons is in agreement 
with the boost invariant Bjorken flow which was done in pUj EH] . 
Only recently we came across a similar work j22] that finds the 
entropy for the apparent horizon. They used a slightly different 
approach than what we did to define the entropy current for the 
apparent horizon. However, it is in complete agreement with our 
results. 

Moreover, a recent work [23j suggests that the apparent horizon 
is a better candidate to describe the dual fluid flow at the boundary. 
This work examines the exact solution of the conformal soliton flow, 
which results in the divergence of the area of the event horizon, while 
the apparent horizon stays finite and constant. This suggests that 
the corresponding entropy of the fluid should be described by the 
apparent horizon. 
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Appendix A 



Notation 



We refer to lower Roman indices that start from the beginning of the 
alphabet a, 6. . . , /i as the whole d + 1 space-time coordinates (bulk 
coordinate). The Greek indices /i, z/... refer to the d dimensions 
hyper-surface of the fluid that "lives" on the boundary of the bulk, 
which are just the indices of the CFT. The lower Roman indices 
that start from i,j... refer to the d — 1 spatial directions of the fluid. 
We commonly use different metrics, which we will specify here: 

• Qab - is the d + 1-dimensional bulk metric. We use this metric 
to lower and raise tensors with bulk indices that are deflned in 
the whole space-time. 

hab - is a (i+ 1-dimensional metric that is restricted to a certain 
hyper-surface. We use this metric to lower and raise tensors 
with bulk indices that are deflned on a certain hyper-surface. 



• 



h^i, - is a (i-dimensional metric of a hyper-surface. We get this 
metric by a projection of Qab into the hyper-surface. We use 
this metric to lower and raise tensors with Greek indices that 
are deflned on a certain hyper- surface. 

Qfj^u - is a (i-dimensional boundary metric. In our calculation 
Qfj^u is the flat metric rjfj,^. Only in the section that deals with 
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the Weyl formalism, do we let it be also a curved metric. 

• rifj,^, - is the rf-dimensional Minkowski boundary metric, with the 
signature convention of "more plus signs" (— , +, +, ...). This is 
the metric of the fluid in the boundary and we will lower and 
raise tensors with boundary Greek indices with that metric, 
i.e., u^ = rj^^Uiy. 

• Qab - is a (i+1-dimensional metric, which defines the d—1 space- 
like hyper-surfaces that create the d-dimensional tube. This 
metric projects the bulk tensor onto the d — 1 hyper- surface. 
We use this metric to lower and raise tensors that we want to 
restrict into the d—1 hyper-surface and that are written in the 
bulk indices, a,b, .... 

• Qij - is a (i — 1-dimensional metric, which defines the d—1 space- 
like hyper-surfaces that create the d-dimensional tube. We use 
this metric to lower and raise tensors that are restricted to the 
d—1 hyper-surface and have i,j... indices. 

We remind the reader of the basic tensors used commonly in this 
thesis: 

The Riemann curvature tensor defined below for bulk vector field 

A", 

VeVrfA'^ - V,V,A« = Kcd A' (A.l) 

^bcd — ■■- bd,c '- cb,d "r -■- ce'- db '- de'- cb l^-^J 

Additionally we define the Ricci tensor and the Ricci scalar by, 

Rbd = Rlad (A- 3) 

R = g^'Rab (A.4) 
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Appendix B 



Detailed Calculations 



For the hyper-surface $ = r — i?, where R is constant, we have the 



basis vectors written in (5.3). Then the metric on the surface, to 
any order, is: 

h^u = gabe^^el = g^,^ (B.l) 

The inverse metrics to the zero and the first order are: 



The Christoffel symbol, with respect to the metric of the hyper- 
surface <l> to the first order, is: 
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p<7(l) 



1 Lo-A 



h"^ {h^\,u + hy\^^ — h^u^x) 



2 {bRY 



U(nd,Au'^ + dunUuP d\ In 6 



H/it/zy) 



^■Uj.r' CA . 



-P^^d. iu^u,) + 



^A i"/i"iy; 



JWf (^^^^"'^^ 



— -dui^^dy) In 6 — du^Ufj_Uud\ In 6 



+m^(9a (m^m,.) 



(B.4) 



fM(2) ^ ^ L 

^"^ d-lRf {bR) '^' 



{dydJnf{bR)-d,dX) 



(B.5) 



T;^Muo 



2Rf{bR) 



-2uPdyd,ux - 2^5, In / {bR) 
+ {l-f{bR))u^{dyd,u^ 
-u^dpU^u^'d^ux - dpU^dxu^) ) (B.6) 



We will specify the detailed calculation of the Christoffel symbols 
for the entire bulk space-time. 

The Christoffel symbols that we can calculate in full to all orders, 
of the full bulk metric gab, are: 

rr 7^9 \9ra,r ~r ()ra,r 9rr,a) 



r)9 \9rr,r ~r Qrr,r 9rr,r) 

~^~^9^ {9r\,r + 9r\,r — 9rr,\) = 



(B.7) 
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^ rr r)9 \9ra,r ~r ()ra,r 9rr,a) 

7)9 \9rr,r ~r ()rr,r 9rr,r ) 
+ 2^''^ i9rX,r + 9r\r " 9rr,\) = (B.8) 



The Christoffel symbols that are needed to calculate the Ricci 
asor, wi 
order are: 



tensor, with respect to the metric g^^ = g^^ + g^^ , to the first 



p!^(0+l) ^ ^ ua ( , _ X 

= 19"'' {9,r,. + 9.r,, - 9,.,r) 



-bdr {r^F {br)) u^cr^^ - -j—u''u^u„ 

1 1 

+ -u''u''dp {upU„) - —dr {r^f (br)) u^UaU^dpu" 

+u'dXP^a (B.9) 



r;:r^^ = l9''''i9,a,r + 9ra,,-9,r,a) 

1 1 

'7^9 \9t^r,r + 9rr,iM 9lJ.r,r) + '^9 {9fiX,r H~ 9r\,^i 9lJ.r,\) 

= lP;: + ^2+^l^ + bdriFibr))a; (B.IO) 
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'- Ij,!/ nil \yiJ,a,u ^ yua,ii ytMi/,a) 

1 1 

= ^r'^f ipr) {-2d(pU^) - dr {-r'^f (br) u^u^ + r'^Pf.u)) 

+ V[29(.(/(6r)..,)) 

-2(9(^M^) - M^(9p (-/ (6r) u^Uy + P^y 



+r^f (br) I - bdr [r'^F {br)) a^^ 
dxu^ 1 



Uf,U,y + -u ux iu^u^) 



d-l^^ 2 



l^rS. (rV (6r)) u,u. + 2r'^P,. (B.ll) 



1 

pr(0+l) ^ :!: m/ , _ N 

-*- /ir r)i/ \yiJ,a,r ^ yra,^i y^ir,a) 

qS* \9l-ir,r ~r Qrr,^ 9fJ.r,r) 

+ lfX^^^'H9,X,r + 9rX,, - 9,r,x) 

= Idr {r'f{br))u,-^u, (B.12) 

The Christoffel symbols that are needed to calculate the expansion 
parameter 6'(£), with respect to the metric 5^^^ = g^f^ +9ab ~^9ab 
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to the second order, are: 

^ liv nil \yiJ,a,u ^ yua,ijL yfMi/,a) 

= ^r'^f ipr) {-2d(pU^) - dr {-r'^f (br) u^u^ + r'^P^y)) 



+r^f (br) ( - hdr [r^F (br)) a 



dxu^ 1 



d- 



-U^Uy + ^u^ux {u^Uy) 



■^rdr {r'f M) u,uy + Sr^l^P,. (B.13) 
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■nr(0+l+2) 
rfj. 



r\9 ydfia^r ~r Qra^fj, 9^ir,a) 

1 1 

'7)9 \9t^r,r + 9rr,iM fl'/ir,rj + 'p:9 \9lJ.X,r + 9r\,ii 9f^r,X) 



dxu^ . 

-u„- — u.,a 



d-1 



fi'-'r 



2{br) 



b i^' i' - "^'0 



a/3 



-2^./3C^"^-(H'^2(6r)"|^ 



+ -u^dxu^P;:dxu^ - brdrF (br) u^d^u^a; 



2r 



+ 



+ 



1 1 



2L (br) 
r {brf-^ 2'^'' \{hrf-^ ) ' r d - 2 



Ig (2Libr) 



xP^{d,-{d-l)A,)a> 



1 
'2r 



u^dxu^d^u. 



(B.14) 



pKO+l+2) ^ i i^a/ , _ ^ 

ofi* \9ur,r ~r yrr,u 9i'r,r) ' ~^9 \9y\,'r ~r QrX^u 9ur,X) 



-{(d-^y'dr^. 



{d— 1) — b'^dr-Ki (br) aapc 



a/3 



+ -P''^u''dpa^pdr {H2 {br) - H, {br}) 
1 



+ —Ual3UJ' 



al3 



(B.15) 
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■p 1/(0+1+2) 



n9 \9va,^ i 9^a,u 9v^,a) 

1 1 



1 



(C5-i)'^''9^(-2m(^V,) + 0.,) 



u„- — - +u„a 



'd-1 



fi'-'r 



b ('■' i' - C'-'l 



2{br 



al3 
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u^dxu^P^^dxu'' + brdrF (br) u^dxu^al 



2r 
+ -^P'idx-id-3)Ax)io 



+ 



rd-2 ^ 

2L{br) 1 l2L{br) 



+ 



r (6r 
1 1 



,d-2 



-9,. 



\d-2 



rrf-2 ' '^ 

\ 

+—u^dxu''d^Uf, 
2r 

_pr{0+l+2) 



2 \ {br 
P^{dx-{d-l)Ax)a 



(B.16) 



^^^.(0+1+2)^^ = 



(B.17) 
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Appendix C 

The Constarint Equations to The Third 
Order 



We will present here a different method than what we presented in 
chapter 6 to find the constraint equations to the third order of the 
event and apparent horizons. 



In order to do so, we look at equation (4.23) and use (4.20) to 



get the following form of the constraint equations: 

Rabm'^e''^ = (C.l) 

In order to check that we get the same constraint equations from 
the two hyper-surfaces, we calculate the difference between the con- 
straint equations evaluated on the event horizon and the constraint 
equations evaluated on the apparent horizon, to the third order in 
derivative expansion, i.e.: 



EH -AH = i?afem"et| - Rabm''e''J (C.2) 

The quantities will be calculated from the full second order bulk 



metric (3.7). Bearing in mind that the difference in m"' only starts 
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from the second order, and that the difference in Rab and e" only 



starts from the third order, we can write (C.2) in the following 
manner: 



EH -AH 









+ m' 



.«{2)| 



I r=rEH 



+ e' 



3^(3) I 



m 

X3) 



a{2)\ 






"an 



+ I m* 
+ (rri 



«{3)| _ rr)''(^)| 



i?i?m^(o) 
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r=rEH 
«{2)| _ ™a{2) 



+ 1 rT 



I r=rEH 



r=rEH 



R 



m 

(2) 






?(0)^b(l) 



I r=rAH 






ab 



r—r-AH 



^a(0)g6(l) 






a(l) 



(C.3) 



where the superscript refers to the order that we are evaluating the 
specific quantities, i.e., if it is written ^^' only the first derivative 



terms should be taken (no zero terms). The quantities in (C.3) that 
are outside the brackets should be evaluated at the location r = b^^, 
for instance: 



m 



r(0)| 



=6-1 



9 



rr{0) I 



=6-1 



r^fibr)\ 



:6-l 



(C.4) 



Here we write some of the quantities that will help us to evaluate 



(C.3) 



Ri'J = 0, 

m"(o) = 0, m"^^^ = 0, 

^M2)| =^M2)| 

\r=rEH \r=rAH ' 

-/~L„,. Hill/. -^ri/ -^Tiy -^n/ ^ "i 



A^) 



rW 



SO) 



^^lE'H/AH ~ ^ime'h/AH ~ ^tJ'''^EH/AH ~ ^M" ' 



(C.5) 



.«{3) 



ef ) = 9,r(2) 
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We know that only the third order terms that will be different 
are the one that are composed of the second order horizon location 
r^'^\ which is different in both hyper-surfaces. For instance, the only 
part that will contribute in the third order Ricci tensor, is the Ricci 
tensor calculated to the first order in derivatives evaluated on one of 
the hyper-surfaces and from this we will take the third order terms 



that come from r 



(2) 
EH/AH- 



We denote = to represent the only terms 



that will not cancel out in the difference m""^^' \ 
for example we will write. 



I r=rEH 



- m"(3) I 



\r=rAH 



m 



/^(3)| 



I r=rEH/AH 



friS) _ g,mQ^ (^,g)^^^^ (^)) _ ^mK2)5^ (,(0) (^)) 



-U 



'Oxuf- 



,P^'^d„r 



(2) 



r 

1 ^. .. 1 



t^' EH/AH 



-P^"'d,.b~^ 






u^dxu" + —b-^u^dxu 



U2pt^iya (2) 

" ^ '^f^' EH/AH 



-b'r'^^H/AH^'dxU^' - ^'P''d,rf^,AH 



(C.6) 



m 



r(3)l 



^'Teh/ah 



fr(3) _ grmg^ ^,m^^^^ (^)) _ gr.i.)Q^ (^(0) (^)) 



O ^^^^ ^^ (2) 

-Ir- : U OyT 



d-1 



A o (2) 

— —2r 



dxu^ 

EH/AHjZrj 



EH/AH 



EH/AH 



(C.7) 



where r eh/ah is the event or apparent horizon full location up to 
the second order, and r*^°) = h^^ . 



Now we are ready to evaluate each term in (C.3): 
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(R^all -421 )m'^W = 

= (4? I -421 )^'^'^ + (r^I\ -421 

y 'A* \r=rEH ^ '"''=''' ah J \ '^f^ ir=rEH i^f^ \r=rAH 

= 



m'^(i) 






r=rEH 






^a(O) Ml) 



(421 -421 )m^^'^ef^ + (R^^\ -421 ) 

\^ rr \r=rEH ^^ ^r=rAH J ^ \ A"^ lr=rB^ i''^ \i~tah } 



m^^'^ef^ 



m'^(2)| -m'^^^)] 42e?'^ 



m'-(2) I - m'-(2) I ) 40) 6^(1) + f m^^^^) I _ ^.(2) I ^ ^(0)gr{l) 



a(3)| _^-(3)| ')i?(0) = 

i'-(3)| -m'-(3)| ) 40)^/^.(3)1 _^K3)| \m 

= -cifr(2)| -r(2)| ]2^u, 

\ \r=rEH \r=rAH J d — 1 

-du^u, (9. (r(^) (x)) l^^^^ - d. (r(^) (x)) |_^ J 
+rfP: f 9, (r(2) (x)) I - d, (r(^) (x)) I ) 

M y V ^ ' / \r=rEH ^ '\i'=fAH) 

= -d f r(2) I - r(2) I ) 1^21^..^ - u'dxu, 

V l'-='-™ \r=rAHj y d-1 '^ 

+4a,(r(^)(x))L_ -9,(r(^)(a:))| 
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UT I - ef ^ I ) 4?m'^°^ + UT I - ^T I ) ^^V'^°^ 

y P \r=TEH P ^f='''AHJ \ ^f=1'EH ^ ^ '''='<' AH J 

-d(a,p)(a:))|_,^^-9,(r(^)(x))|_J 



f^-(2)| _^a(2)| A^(l)^ 

V \r=TEH \r=TAHj °-f^ 

= f m^(2) I _ ^.(2) I \ ^(1) /^r.(2) I _ ^K2) I ^ ^(1) ^ 

V ir=rEH ir=rAHj '^f^ \ ir=rEH ir=rAH J '^'^ 



(421 -<1 )m'^(°) = 

\^ rfJ. \r=rEH '^ ''''=''' AH J \ '^f^ 'r=rEH ^ 

\^ '^fJ.lr=rEH '^fJ-lr=rAH/ 

= d f r(2) I - r(2) I ) (2^u, - u'^^u^ 



i^M I r=rAH , 



We seek the third order constraint equations. This is the reason 
that in our calculation we used the previous constraint equations 



(5.16), (5.17) that are of the second order. Assembling all of the 



above terms and inserting them into (C.3), reveals that our state- 
ment is correct, for we found: 



EH -AH = /2a6mV| - i^ab^V I =0 (C.8) 

ao tJ.\r=rEH IJ-\r=rAH ^ ' 
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